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The problem of brain connectivity has 
been gaining more and more interest in 
the last years. Connectivity can be esti- 
mated by different techniques and at the 
different levels of the hierarchy of the ner- 
vous system. Here we shall consider func- 
tional connectivity at the level of the brain 
structures, in particular connectivity mea- 
sures derived from scalp EEC measure- 
ments. A multitude of estimators are being 
used for connectivity estimation: linear, 
nonlinear, bivariate and multivariate. The 
disadvantage of bivariate estimators is con- 
nected with the fact, that in virtue of com- 
mon feeding many spurious connections 
are found. This fact was demonstrated in 
works: (Blinowska et al., 2004; Kus et al., 
2004). Let us consider the common sit- 
uation when a source is emitting activity 
measured at N electrodes. In case of bivari- 
ate measures beside N true connections 
also false connections between all elec- 
trodes will be find in virtue of common 
feeding and their number will be [N(N — 
l)/2 — N], so they will outnumber the 
true connections. 

This fact directed the interest of the 
scientific community to the multivariate 
methods of connectivity estimation. In 
particular, multivariate estimators based 
on the Granger causality principle allow 
to estimate directed connectivity and 
does not produce spurious connections. 
Among these estimators widely applied 
are: Directed Transfer Function (Kaminski 
and Blinowska, 1991) and Partial Directed 
Coherence (Baccala and Sameshima, 
2001) based on the Multivariate 
Autoregressive Model (MVAR). Currently 
DTF and PDC are commonly used and 



have became a part of various signal 
processing packages, e.g., eConnec- 
tome (http://econnectome.umn.edu) , 
Octave-Forge (http://octave.sourceforge. 
net/tsa function mvfreqz.html) , Epilab 
(http://www.epilepsiae.eu/project_outputs/ 
epilab _ software), SIFT (http://sccn.ucsd. 
edu/wiki/SIFT). 

Unfortunately we have noticed that 
quite often the application to DTF of inap- 
propriate preprocessing routines produces 
misleading results. Moreover, the use of 
these routines is not necessary in virtue 
of the properties of DTF. As the authors 
of the DTF estimator, we are particularly 
concerned that the method is applied cor- 
rectly. 

Among the pre-processing methods 
used before DTF application the most 
common approach is to project the signals 
into the source space in order to eliminate 
volume conduction effect. However, DTF 
for a pair of channels i and j is nonzero 
only, if there is a phase difference between 
channels i and j. 

The volume conduction is a propaga- 
tion of the electromagnetic field, so it does 
not produce a phase difference at elec- 
trodes. Therefore, DTF is practically insen- 
sitive to volume conduction. The fact that 
the estimators of connectivity based on the 
phase difference between channels are not 
influenced by the volume conduction was 
also recognized in Stam et al. (2009). The 
influence of a signal of a constant phase on 
DTF was demonstrated by means of simu- 
lation. For the set of EEC signals we have 
added a sinusoid, with the same phase for 
each signal. The amplitude of that sinusoid 
(of 20 Hz) was similar to the amplitudes 



of EEC signals. The result is shown in 
Figure 1. We can observe a prominent 
peak at 20 Hz in the power spectra of the 
signals, but this peak is absent in DTF 
functions. 

The fact that DTF is not influenced 
by volume conduction effects is further 
supported by the excellent topograph- 
ical agreement of DTF results with 
the evidence obtained from anatomical 
data, imaging studies and physiologi- 
cal experiments. Clear cut patterns of 
propagation emerged from these studies. 
As the examples may serve the stud- 
ies of finger movement (Ginter et al., 
2001; Kus et al, 2006), localization of 
epileptic focus (Franaszczuk et al., 1994) 
results obtained for Constant Attention 
Test (Blinowska et al., 2010) and work- 
ing memory paradigm (Brzezicka et al., 
2010; Blinowska et al, 2013). More 
examples may be found in Blinowska 
(2011); Blinowska and Zygierewicz (2011). 
The animations illustrating dynamically 
changing patterns of connectivity obtained 
by the Short-time Directed Transfer 
Function (SDTF) in some of the above 
mentioned experiments are available at 
http://brain.fuw.edu.pl/~kjbli. 

Since DTF is immune to volume con- 
duction the pre-processing procedures 
such as projections on the cortex sur- 
face or Laplace transform are not needed. 
Moreover, they destroy the original corre- 
lation structure of the set of signals. If we 
mix the information from different chan- 
nels calculating for example Laplacian, 
we influence the correlation between sig- 
nals and the information about the phase 
relations between channels is disturbed. 
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FIGURE 1 | The DTFs (as functions of frequency) estimated for a set of 
EEG signals. (A) The results for original dataset, (B) same data with 20 Hz 
sinusoid added with constant phase to each channel. On the diagonals of the 



panels power spectra are shown. The propagation from the channel marked 
below the pictures to the channel marked at the left. At the (C) EEG from one 
channel with superimposed sinusoid. 



As a result, information on causality or, 
in another words, on the propagation of 
activity from channel; to i is lost. 

The results of the works where the 
preprocessing involving projection on 



the cortex was applied show a disorga- 
nized structure of connectivity. The pre- 
processing of the data before the DTF 
application should be limited to subtrac- 
tion of the mean, possibly division by 



the variance and digital filtering. However, 
the filtering must not influence phases of 
signals. These can be achieved by filter- 
ing forward and backward (e.g., Matlab 
procedure filtfilt). The signals should be 
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referenced to the "neutral" derivation, e.g., 
linked ears, nose, or similar. No com- 
mon average or bipolar reference may 
be used. There should be no preprocess- 
ing by means Hjorth or Laplace trans- 
form, Independent Components Analysis 
or projection of signals on the cortex. In 
the view of robustness of DTF in respect to 
volume conduction, these procedures are 
obsolete, furthermore, they are harmful 
and may produce misleading results. 

It is worth mentioning that in gen- 
eral methods of connectivity estimation 
based on phase differences are insensitive 
to volume conduction. For instance this 
property holds also for Partial Directed 
Coherence, but not for ordinary coher- 
ence, because it contains common compo- 
nent, which includes activity of no phase 
difference between electrodes. In conclu- 
sion we would like to underline that the 
preprocessing of the signals for estimation 
of connectivity should be based on the 
thorough understanding of the properties 
of the applied methods. 
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